GRB 000418: A Hidden Jet Revealed? 
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ABSTRACT 

We report on optical, near-infrared and centimeter radio observations of 
GRB 000418 which allow us to follow the evolution of the afterglow from 2 to 
200 days after the 7-ray burst. In modeling these broad-band data, we find that 
an isotropic explosion in a constant density medium is unable to simultaneously 
fit both the radio and optical data. However, a jet-like outflow with an opening 
angle of 10-20° provides a good description of the data. The evidence in favor 
of a jet interpretation is based on the behavior of the radio light curves, since 
the expected jet break is masked at optical wavelengths by the light of the host 
galaxy. We also find evidence for extinction, presumably arising from within the 
host galaxy, with Ay os *=0.4 mag, and host flux densities of Fr = 1.1 fiJj and 
F K = 1.7 /xJy. These values supercede previous work on this burst due to the 
availability of a broad-band data set allowing a global fitting approach. A model 
in which the GRB explodes into a wind-stratified circumburst medium cannot 
be ruled out by these data. However, in examining a sample of other bursts (e.g. 
GRB 990510, GRB000301C) we favor the jet interpretation for GRB 000418. 
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1. Introduction 

GRB 000418 was detected on April 18, 2000, at 09:53:10 UT by the Ulysses, KONUS- 
Wind and NEAR spacecraft, which are part of the third interplanetary network (IPN) 
. The event lasted ~30 s, and a re-analysis of the early Ulysses data ( |Hurley, Cline &j 



Mazets 2000]) gives a fluence of 4.7 x 10~ 6 erg cm" 2 in the 25-100 keV band. A fit to the 
total photon spectrum from the KONUS data 2 x 10~ 5 erg cm -2 in the energy range 15 
- 1000 keV. Intersecting IPN annuli resulted in a 35 arcmin 2 error box, in which Klose et 
al. (2000b) identified a variable near-infrared (NIR) source. The early R-band light curve 
of this source was described by Mirabal et al. (2000) as having a power-law decay t~ ' 84 , 
typical for optical afterglows. The redshift for the host galaxy of z ~ 1.119 was measured 
by Bloom et al. (2000) from an [Oil] emission line doublet. Assuming cosmological 
parameters of Qm=0.3, A =0.7 and H =65 km s _1 Mpc -1 , this redshift corresponds to a 
luminosity distance dL = 2.5 x 10 28 cm and gives an implied isotropic 7-ray energy release 
of E 7 = 1.7 x 10 52 erg. 

Klose et al. (2000b) have recently summarized optical/NIR data observations of 
GRB 000418. In this paper we present additional optical/NIR data and a complete set of 
radio observations between 1.4 GHz and 22 GHz, from 10 to 200 days after the burst. We 
use this broad band data set to fit several models, deriving the physical parameters of the 
system. 



2. Observations 

2.1. Optical Observations 

In Table |1] we present deep optical photometry obtained at Palomar, KeckP], and 
MDM observatories covering six weeks following the GRB as well as data from the extant 



12 The W. M. Keck Observatory is operated by the California Association for Research in Astronomy, a 
scientific partnership among California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration. 
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literature. 

All of the optical data was overscan corrected, flat-fielded, and combined in the usual 
manner using IRAF ( [Tody 19931 ). PSF-fitting photometry was performed relative to several 
local comparison stars measured by Henden (2000) using DoPhot ( (Schechter, Mateo jg 
[5aha 1993|) . Short exposures of the field in each band were used to transfer the photometry 
( |Hcnden 2000|) to several fainter stars in the field. Several of the Keck+ESI measurements, 



and the Palomar 200" measurement were made in Gunn-r and Gunn-i respectively and were 
calibrated by transforming the local comparison stars to the Gunn system using standard 
transformations ( |Wade et al. 1979^ , [Jorgensen 1994] ). We add an additional 5% uncertainty 



in quadrature with the statistical uncertainties to reflected the inherent imprecision in these 
t r ansf or mat ions . 

The Ks-band image of the field was obtained on the Keck I Telescope on Manua Kea, 
Hawaii with the Near Infrared Camera. We obtained a total of 63 one-minute exposures 
which we reduced and combined with the IRAF/DIMSUM package modified by D. Kaplan. 
There was significant cloud and cirrus cover and so the night was not photometric. 

The HST STIS/Clear image was obtained on 4 June 2000 UT as part of the TOO 
program # 8189 (P.I. A. Fruchter) and made public on 2 September 2000 UT. Five images 
of 500 s each were obtained which we combined using the IRAF/DITHER task. The final 
plate scale is 25 milliarcsec pixel -1 . 

We corrected all optical measurements in Table |I] for a Galactic foreground reddening 
of E(B - V) = 0.032 ( [Schlegel, Finkbeiner fc Davis 199g| ) at the position of the burst 



(/, b) = (261.16,80.78) before converting to flux units ( |Fukugita, Shimasaku fc Ichikawa 
[19951 [Bessell fc Brett 1988Q assuming R v =3.1. 



2.2. Radio Observations 

Radio observations were undertaken at a frequency of 15 GHz with the Ryle Telescope. 
All other frequencies were observed with either the NRAO0 Very Large Array (VLA) or 
the Very Long Baseline Array (VLB A). A log of these observations can be found in Table |[ 
The data acquisition and calibration for the Ryle and the VLA were straightforward (see 
Frail et al. 2000a for details). 



13 The NRAO is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. NRAO operates the VLA and the VLBA 
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The single VLBA observation was carried out at 8.35 GHz with a total bandwidth 
of 64 MHz in a single polarization using 2 bit sampling for additional sensitivity. The 
nearby (<1.3°) calibrator J1224+2122 was observed every 3 minutes for delay, rate and 
phase calibration. Amplitude calibration was obtained by measurements of the system 
temperature in the standard way. The coordinates for GRB 000418 derived from the VLBA 
detection are (epoch J2000) a = 12 /l 25 m 19.2840 s (±0.015 s ) 5 = +20°06'11.141 // (±0.001"). 



3. The Optical Light Curve and Host Galaxy 

In Figure [I] we display the R and K-band light curves constructed from measurements 
in Table |l[ The pronounced flattening of the R-band light curve at late times is reasonably 
attributed to the optical afterglow fading below the brightness of the underlying host galaxy. 
A noise-weighted least squares fit was made to the data of the form f R = f t% + fhost for 
which we derive f a = 23.4 ± 2.1 jiJy, a Q = -1.41 ± 0.08 and f host = 1.08 ± 0.06 fiJy with a 
reduced Xr — 0.94. Our inferred R-band magnitude for the host galaxy Rhost = 23.66 ± 0.06 
is nearly identical to that obtained from a similar analysis by Klose et al. (2000b). In order 
to estimate the effect of the host in other optical bands we scaled Rhost for GRB 000418 
to a spectrum of the host galaxy of GRB 980703 ( jBloom et al. 199"E| ) (z = 0.966) whose 



magnitude was measured in seven broad-band colors (B, V, R, I, J, H, and K). Our results 
indicate that 50-100% of the flux in some bands is due to the host galaxy after the first 
10 days. Therefore, for the afterglow modeling in §|5] we chose not to include the late-time 
measurements of GRB 000418 in the B, V, and Gunn-i bands. 



4. The Radio Light Curves 

In Figure [l] we display the radio light curves at 4.86, 8.46, 15 and 22 GHz. To first 
order all four frequencies show a maximum near 1 mJy on a time scale of 10 to 20 days. 
There is no discernible rising trend at any frequency. This is most clear at 8.46 GHz, where 
beginning 10 days after the burst, the light curve undergoes a steady decline, fading from 1 
mJy to 0.1 mJy over a six month period. The temporal slope of the 8.46 GHz light curve 
after the first two months a m d = —1.37 ± 0.10 (xl = 1-4) is similar to the optical R-band 
curve a opt = —1.41 ± 0.08. 

Superimposed on this secular decrease, there exist point-to-point variability of order 
50%, especially in the early measurements. We attribute these variations to interstellar 
scintillation (ISS) QGoodman 1997 , [Walker 1998|) . The method by which we estimate the 
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magnitude of the intensity fluctuations induced by ISS as a function of frequency and time 
is described in full by Berger et al. (2000). Briefly, we estimate the magnitude of scattering 
with the model of Taylor & Cordes (1993), and use this to calculate the transition frequency 
vq between weak and strong scattering using Walker (1997). The normalizations used in 
Goodman (1997) give slightly larger values of u . 

In the direction toward GRB 000418 we derive uq ~ 3.6 GHz and therefore most of 
our measurements were taken in the weak ISS regime. In this case the modulation scales 
as z/~ 17 / 12 , with a maximum of 65% expected at 4.86 GHz and 30% at 8.46 GHz. At 15 
GHz and 22 GHz we estimate that the ISS-induced fluctuations are only a fraction of 
the instrumental noise. The expansion of the fireball will eventually quench ISS when 
the angular size of the fireball exceeds the angular size of the first Fresnel zone at the 
distance of the scattering screen. The fireball size, and hence the quenching timescale, is 
model-dependent, and we use the derived total energy and density from the global fits 
(see §|5| below) to estimate this time for each model. For example, in a simple spherical 
fireball this occurs after 15 days at 4.86 GHz and 10 days at 8.46 GHz, and thereafter the 
modulation indices decline as t -35 / 48 . We note that the observed fluctuations at 4.86 and 
8.46 GHz conform to the predicted level of ISS, but that the measurements at 8.46 GHz 
from around 50 days after the burst deviate by a factor of three from the predicted ISS 
level. 

In addition, we use the scintillation pattern to estimate the true xl f° r eacn model, 
by adding in quadrature to the instrumental noise an additional ISS-induced uncertainty, 
ciss = mpFv,iDodeh where m p and -F^modei are the modulation index and model flux density 
at frequency v, respectively (Berger et al. 2000). 



5. Global Model Fits 

The optical and radio data presented here have allowed us to track the evolution of the 
GRB 000418 afterglow from 2 to 200 days after the burst. With careful modeling of the light 
curves, it should be possible to infer the physical parameters of the blast wave and thereby 
gain some insight into the nature of GRB progenitors. In particular, the hydrodynamic 
evolution of the shock is governed by the energy of the explosion, the geometry of the 
expanding ejecta shock and the type of environment into which the GRB explodes ( |SariJ 
piran fc Narayan 1998| , |Wijers fc Galama 1999| , |Chevalier fc Li 1999j [Panaitescu fc Kumar| 



2000a| ). We will consider three basic models: a spherical explosion in a constant density 



medium, collimated ejecta (i.e., jets) interacting with a constant density medium, and a 
spherical explosion in a wind-blown medium. 
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The starting point for any afterglow interpretation is the cosmological fireball model 
(e.g. [Meszaros fc Rees 1997| , |Waxman 1997| ). A point explosion of energy E expands 



relativistically into the surrounding medium (with density p oc r~ s , where s = for constant 
density ISM and s = 2 for a wind) and the shock produced as a result of this interaction is 
a site for particle acceleration. The distribution of electrons is assumed to be a power-law 
of index p, and the fraction of the shock energy available for the electrons and the magnetic 
field is e e and e#, respectively. The values of these three quantities (p, e e and e#) are 
determined by the physics of the shock and the process of particle acceleration and in the 
absence of detailed understanding are taken to be constant with time. 

The instantaneous broad-band synchrotron spectrum can be uniquely specified by the 
three characteristic frequencies v a , v m , and v c (i.e., synchrotron self-absorption, synchrotron 
peak, and cooling), the peak flux density f m , and p. For this work we adopt the smooth 
spectral shape as given by Granot et al. (1999a, 1999b) rather than the piecewise, broken 
power-law spectrum used by other authors (e.g., |Wijers fc Galama 1999Q . The evolution 
of the spectrum (and thus the time dependence of v a , v m , v c and f m ) is governed by 
the geometry of the explosion (spherical or a collimated into a jet-like outflow), and the 
properties of the external environment (constant density or a radial density profile). Our 
approach is to adopt a model (sphere, wind, jet, etc.) and solve for the above spectral 
parameters using the entire optical and radio data set. The advantages and details of global 
model fitting are discussed by Berger et al. (2000). 

The simplest model is a spherically symmetric explosion in a constant density 
medium (ISM: Sari et al. 1998). The total xl f° r this model (see Table gives a highly 
unsatisfactory fit to the data. On close inspection (Figure P we find that the model 
systematically underpredicts the optical flux. Adding extinction from the host galaxy 
only makes this worse. The fundamental difficulty with the ISM model is that it predicts 
f m = constant, independent of frequency. In this case, since it is the radio data that is 
responsible for defining the peak of spectrum, it results in a value of f m that is too low at 
higher frequencies. 

To obtain better fits to the joint optical and radio data sets we look to models for 
which f m is time- dependent. One such model is a collimated outflow into a medium with 
uniform density ( Jet: |Rhoads 1997| , [Rhoads 1999| , |Sari, Piran fc Halpern 1999| ) . The clearest 



observational signature of the Jet model is an achromatic break in the light curves at tj (e.g. 
Harrison et al. 1999|). At radio wavelengths (i.e., below u m ) at tj we expect a transition 



from a rising; t 1/2 light curve to a shallow decay of t 1/1,3 , while at optical wavelengths the 
decay is expected to steepen to t~ p . These decay indices refer to the asymptotic values. 



Detecting a jet transition at optical wavelengths may be difficult if it occurs on 
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timescales of a week or more. In these cases the afterglow is weak and the light from the 
host galaxy may start to dominate the light curve (e.g., Halpern et al. 2000). In such 
instances radio observations may be required to clarify matters, since the radio flux is 
increasing prior to tj and changes in the lightcurve evolution due to the jet break are easily 
detected. Indeed, the jet in GRB 970508, which was very well observed in the radio is not 
discernible in the optical data. In this case, Frail, Waxman & Kulkarni (2000) found a 
wide-angle jet with an opening angle of 30° and tj ~ 30 days (but see |Chevalier fc Li 20001 ). 

A Jet model with tj p» 26 days fits the data remarkably well (see Figure [I]). The 
strongest point in favor of the Jet model is that it reproduces the broad maximum (~1 
mJy) seen from 5 GHz to 22 GHz. We expect such a plateau at tj as all light curves for 
v a < v < v m reach their peak fluxes (with only a weak v 1 ^ frequency dependence) before 
undergoing a slow decline. Most other models predict a strong frequency dependence in 
peak flux which is not seen in this case. 

Knowing tj and the density of the ambient medium n from the model fit (Table |^) 
we can make a geometric correction to the total isotropic energy E 7 , as determined from 
either the observed 7-ray fluence or the total energy of the afterglow E52, from the fit to the 
afterglow data (Sari et al. 1999, Livio fc Waxman 2000 ). This approach gives values for the 



jet opening angle 9j between 10° and 20°, which for a two-sided jet reduces the GRB energy 
to ~ 10 51 erg. The rapid lateral expansion of the jet also accelerates the transition to the 
non-relativistic regime, resulting in a change in the evolution of the light curves. Since this 
occurs on a timescale £nr ~ tflj 2 ~ 350 days ( [Waxman, Kulkarni fc Frail 1998|) , we do not 



expect the non-relativistic transition to be important for our data. 

There is some freedom in our choice of v c . We know that a cooling break (i.e., 
Aa = —0.25) is not apparent in the R band light curve on timescales of 2-10 days, so we 
searched for solutions with v c above or below this frequency. We found that physically 
consistent solutions (i.e., with non-negative host fluxes, and €3 < 1) were only possible for 
values of v c below the optical band. 

As part of the fitting process we also solved for the host flux density in the R and K 
bands and for any local dust obscuration, assuming an LMC-like extinction law. This yields 
/host (R)= 1.1 /xJy, f host (K)=1.7 /iJy and Ay ost =0A (in the host galaxy restframe). Klose 
et al. (2000b) argued for significant dust extinction with A v osi =0.96. However, they likely 
overestimated Ay ost since they assumed a spherical fireball model and arbitrarily located v c 
above the optical band. Moreover, we find that there is some covariance between the values 
of Ay 05 * and p so that only with a global fit, in which p is constrained by the radio data as 
well as the optical data, we can solve for Ay ost in a self-consistent manner. 
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In view of the claims linking GRBs with the collapse of massive stars flGalama et al. 



|1998| , [Bloom et al. 199~9| , [Reichart 1999| , |Piro et al. 2000| ), we considered a final model of 
a spherical explosion into a wind-blown circumburst medium ( Wind: |Chevalier fc Li 1999| , 
Li fc Chevalier 1999 ). The Wind model (Figure [T]) fits the data as well as the Jet model. 
In fact, the two models have sufficiently similar \ 2 t° be indistinguishable. The close 
match between the temporal slopes of the late-time 8.46 GHz light curve and the early R 
band light curve (see §|J) is a point in favor of the Wind model since a steeper decline is 
expected for a jet geometry. Our failure to distinguish between these two models can be 
blamed on the absence of radio measurements (particularly at millimeter wavelengths) at 
early times. The rapid rise of the flux density below v a and v m in the Wind model and 
the strong frequency dependence of the peak flux (see Figure [I]), make such measurements 
advantageous ( jPanaitescu fc Kumar 2000aD . Moreover, in principle the Wind model can 
be distinguished from the other models by the fact that in this model v c is increasing 
with time (y c oc t 1 / 2 ). However, in this case since v c lies below the optical/IR bands, this 
behavior would be distinguishable only at late time when the host flux dominates over the 
OT. As before we solved for the host flux and any dust extinction, yielding f host (R)=0.9 
fiJy, fh ost (K)=1.3 nJy and Ay ost =0.3. In our view further (and more sophisticated) model 
fits are not warranted by these data. 



In summary, we find that the radio and optical/NIR observations of the afterglow 
emission from GRB 000418 can be fit by two different models. The close similarity between 
the results of the Wind and Jet models has been noted for other GRBs: GRB 970508 (Frail 
et al. (2000), |Chevalier fc Li 2000|) , GRB 980519 flFrail et al. 2000U IChevalier fc Li 1999|) , 
GRB 000301C QBerger et al. 200q , [Li fc Chevalier 2000|) , and GRB 991208 flGalama et al 



2000| , [Li fc Chevalier 2000| ). The resolution of this conflict is important, since it goes to the 



core of the GRB progenitor issue. If the GRB progenitor is a massive star then there must 
be evidence for a density gradient in the afterglow light curves, reflecting the stellar mass 
loss that occurs throughout the star's lifetime ( |Chevalier fc Li 1999| , |Panaitescu fc Kumar 
2000a| ). At present, an unambiguous case for a GRB afterglow expanding into a wind has 
yet to be found. On the contrary, most afterglows are better fit by a jet expanding in a 
constant density medium (e.g., Harrison et al. 1999| , [Halpern et al. 2000| , [Panaitcscu fc 



Kumar 2000b|) and thus we are faced with a peculiar situation. While there is good evidence 



linking GRBs to the dusty, gas-rich environments favored by hypernova progenitors flBloom 



Kulkarni fc Djorgovski 2000| . C^alama fc Wijers 2000]) , the expected mass loss signature is 
absent (or at best ambiguous) in all afterglows studied to date. 
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Table 1. Optical/NIR Observations of GRB 000418 



TIT Date 


Instr. a 


Band 


Mag. b 


Err 


Ref. c 


Apr 20.89 


TNG 3.5m 


R 


21.54 


0.04 


2 


Apr 20.90 


CA 3.5m 


K' 


17.49 


0.5 


2 


Apr 20.93 


CA 1.2m 


K' 


17.89 


0.2 


2 


Apr 21.15 


MDM 2.4m 


R 


21.66 


0.12 


1 


Apr 21.86 


LO 1.5m 


R 


21.92 


0.14 


2 


Apr 26.316 


USNO 1.3m 


R 


22.65 


0.20 


2,4 


Apr 27.26 


MDM 


R 


22.77 


0.23 


1 


Apr 28.170 


P200 


R 


22.97 


0.06 


1 


Apr 28.3 


MDM 


R 


22.86 


0.09 


1 


Apr 28.413 


Keck/ESI 


R 


23.05 


0.05 


1 


Apr 29.26 


MDM 


R 


22.95 


0.11 


1 


May 2.274 


Keck/ESI 


Gunn-i 


23.38 


0.05 


1 


May 2.28 


MDM 


R 


23.19 


0.12 


1 


May 2.285 


Keck/ESI 


B 


24.31 


0.08 


1 


May 2.31 


USNO 1.3m 


R 


23.11 


0.130 


2 


May 3.26 


USNO 1.3m 


R 


23.41 


0.160 


2 


May 4.44 


UKIRT 3.8m 


K 


20.49 


0.40 


2 


May 6.42 


Keck/LRIS 


R 


23.48 


0.10 


7 


May 8.89 


TNG 


R 


23.30 


0.05 


2 


May 8.92 


TNG 


V 


23.92 


0.07 


2 


May 9.82 


USNO 1.0m 


R 


23.37 


0.21 


2 


May 23.93 


TNG 


R 


23.37 


0.10 


2 


May 29.228 


P200 


R 


23.66 


0.15 


1 


Jim 2.88 


CA 3.5m 


R 


23.32 


0.08 


2 


Jun 2.91 


TNG 


R 


23.57 


0.05 


2 



a CA 3.5m=Calar Alto 3.5-meter, USN01.3m = U.S. Naval 
Observatory Flagstaff Station 1.3-meter, ESI=W.M. Keck 
Observatory Echellette Spectrograph-Imager, LRIS=W.M. Keck 
Observatory Low-Resolution Imaging Spectrograph 

b Optical photometry is on the Kron-Cousins and Gunn 
systems and referred to that of Henden (2000.) Data are 
corrected for Galactic extinction corresponding to E(B — V) = 
0.032 derived from the maps of Schlegel et al. (1998). 



l=this work, 2=Klose et al. 2000, 3=Mirabal et al. 200C 



4= |Hcndcn et al. 200C| , 5= |Mirabal, Halpcrn fc Wagner 2000a , 



6=Mirabal, Halpern & Wagner 20001, 7=Metzger & Fruchter 



200C 
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Table 2. Radio Observations of GRB 000418 



Epoch Telescope 




S±cr 


(UT) 




(GHz) 


OJy) 


9000 A™- 98 7^ 
zuuu i\pi ZO. / 


rtyie 


15.0 


550±600 


9000 Any 90 07 

zuuu Apr zy.u/ 


VT A 


8.46 


856±33 


9000 Ar^r ^0 07 
ZUUU -rt.pl ou.u/ 


VT 4 


8.46 


795±37 


9000 A™- ^0 7^ 
zuuu rtpi ou. / o 




15.0 


1350±480 


9000 A/Ta-tr 1 0£ 

zuuu iviay i.ud 


VT A 


4.86 


110±52 


9000 A/Ta-tr 1 0£ 

zuuu iviay i.uo 


vt A 


8.46 


684±48 


9000 A/Tott 9 OQ 

zuuu iviay z.yo 


T?irl£l 

Kyle 


15.0 


850±300 


9000 A/Ta-tT- *3 0/1 

zuuu iviay o.U4 


vt A 


4.86 


1120±52 


9000 A/Tott- Q 0/1 

zuuu iviay o.u4 


vt A 


8.46 


1240±46 


9000 A/Tott- ^ 0/1 

zuuu iviay o.U4 


VT A 


22.46 1100±150 


9000 Ma-tr /L 07 
ZUUU Ivldy 4.y/ 


vt A 


1.43 


210±180 


9000 A/Tott- A 07 

zuuu iviay 4.y/ 


vt A 


4.86 


710±47 


9000 Matr /L 07 
zuuu ividry 4.y/ 


vt A 


8.46 


1020±53 


9000 A/Tcnr /l 07 

zuuu iviay 4.y* 


vt a 


22.46 


860±141 


9000 A/TcnT- 7 1 ft 

zuuu iviay <.io 


vt r a 

V 1 . 1 > A 


8.35 


625±60 


9000 A/Ta-tT- Q 9c; 

zuuu iviay y.zo 


VT A 


8.46 


926±53 


9000 A/Tqtt 1^1^ 

zuuu iviay 10.10 


VT A 


8.46 


963±34 


9000 A/fi^r 1 R 9/1 
zuuu iviay 10.Z4 


vt A 


4.86 


567±50 


9000 A/Tatr 1 ft 9/1 

zuuu iviay io.z4 


VT A 


8.46 


660±50 


9000 A/Te^r 18 9/1 
zuuu iviay 10.Z4 


vt a 


22.46 


610±114 


9000 A/fa-*?- 99 91 
zuuu iviay zz.zi 


VT A 


8.46 


643±38 


9000 A/fcnr 9fi 09 
zuuu iviay zu.yz 


vt a 


4.86 


1105±51 


9000 A/Ta-tT- 0f\ 09 

zuuu iviay zu.yz 


VT A 


8.46 


341±50 


9000 Inn 11/1 

zuuu jun i.i^: 


vt A 


8.46 


556±43 


9000 Inn 11/1 

zuuu jun 1.14 


vt A 


22.46 


710±16 


9000 Inn *3 0/1 

zuuu jun o.U4 


\H A 


8.46 


517±34 


9000 Tnn 7 01 

zuuu jun ( .ui 


vt A 


8.46 


238±38 


9000 Tnn 110^ 
ZUUU J LL11 J_ J_.czO 


vt A 

V i . . \ 


8.46 


230±33 


2000 Jun 15.13 


VLA 


8.46 


325±30 


2000 Jun 20.10 


VLA 


8.46 


316±30 


2000 Jun 23.19 


VLA 


8.46 


306±29 


2000 Jun 27.08 


VLA 


8.46 


296±22 


2000 Jul 02.98 


VLA 


8.46 


274±22 


2000 Jul 10.04 


VLA 


8.46 


178±24 


2000 Jul 22.81 


VLA 


8.46 


152±23 


2000 Jul 22.81 


VLA 


4.86 


192±25 


2000 Jul 28.50 


VLA 


8.46 


168±22 


2000 Jul 28.50 


VLA 


4.86 


191±25 


2000 Aug 17.74 


VLA 


8.46 


119±25 


2000 Aug 17.74 


VLA 


4.86 


235±31 


2000 Aug 21.65 


VLA 


4.86 


142±35 
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Table 2 — Continued 



Epoch Telescope 


^0 


S±cr 


(UT) 


( 


GHz) 


(/*Jy) 


2000 Aug 21.65 


VLA 


8.46 


87±31 


2000 Aug 25.78 


VLA 


4.86 


238±34 


2000 Aug 25.78 


VLA 


8.46 


166±27 


2000 Aug 27.89 


VLA 


8.46 


100±25 


2000 Sep 10.73 


VLA 


8.46 


148±25 


2000 Sep 18.68 


VLA 


8.46 


55±20 


2000 Sep 26.62 


VLA 


8.46 


85±22 


2000 Nov 06.55 


VLA 


8.46 


94±14 


Note. - The 


columns 


are 


(left to 



right), (1) UT date of the start of 
each observation, (2) telescope name, (3) 
observing frequency, and (4) peak flux 
density at the best fit position of the 
radio transient, with the error given as 
the root mean square noise on the image. 
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Table 3. Model parameters for GRB 000418 



Parameters a 


ISM 


Jet 


Wind 


Va (Hz) 


4.1xl0 9 


1.7xl0 9 


30xl0 9 


v m (Hz) 


2.3xlO n 


1.8xl0 10 


5.8xlO n 


v c (Hz) 


2 x 10 15 


10 14 


1.8xl0 13 


f m (mJy) 


2.5 


3.4 


10.4 


P 


2.3 


2.4 


2.2 


tj e t (days) 




25.7 




A host 


0.0 


0.4 


0.3 


X 2 /dof 
n or A* 


326/54 
11 
0.01 


165/54 
10 
0.02 


184/54 
4 
0.14 



e B 0.05 0.06 0.04 

e e 0.03 0.10 0.07 



a For the ISM and Wind models u a , u m , v c 
and f m are the self-absorption, synchrotron peak, 
and cooling frequencies, and the peak flux density, 
respectively on day 1. For the Jet model these 
values are referenced instead to the jet break time 
tj=25.7 d. p is the electron power-law index and 
A v is the V band extinction in the rest frame 
of the host galaxy (2= 1.1 18), assuming an LMC- 
like extinction curve. Note that in the Jet and 
ISM models v c was fixed at 10 14 and 2 x 10 15 Hz, 
respectively, but in the Wind model it was left as a 
free parameter. The resulting values of \ 2 include 
an estimated contribution of interstellar scattering 
(ISS) and the increased error in subtracting off a 
host galaxy flux from each of the optical points. 
The model parameters arc the total isotropic 
energy E 52 in units of 10 52 erg, the ambient density 
n in cm~ 3 or in the case of the Wind model the 
parameter A* as defined by Chevalier & Li (1999). 
e e and eb are the fraction of the shock energy in 
the electrons and the magnetic field, respectively. 
The true uncertainties in the derived parameters 
are difficult to quantify due to covariance, but we 
estimate that they range from 10 — 20% 
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Fig. 1. — Radio and optical light curves for GRB 000418. The observing frequency (or band) 
is shown in the upper right corner of each panel. Optical magnitudes were first corrected 
for Galactic forground reddening before converting to flux units. For display purposes the 
R band flux densities have been increased by a factor of 10. The 8.46 GHz measurements 
on August 25 and September 18 are 3-epoch averages taken over a period of 7 days and 15 
days, respectively. The dotted, solid and dashed lines are light curves for ISM, jet and wind 
models, respectively. They were derived from a global fit to the entire broad-band dataset. 
See text for more details. 



